Ion size effects on the electric double layer of a spherical particle in a 
realistic salt-free concentrated suspension 



O 
(N 

>> 



Rafael Roa/ Felix Carrique,^'* and Emilio Ruiz-Reina^ 

^ Ftsica Aplicada I, Umversidad de Malaga, Spain, 
^ Fisica Aplicada II, Universidad de Malaga, Spain. 
(Dated: January 20, 2013) 

A new modified Poisson-Boltzmann equation accounting for the finite size of the ions valid for 
realistic salt-free concentrated suspensions has been derived, extending the formalism developed for 
pure salt-free suspensions [Roa et al., Phys. Chem. Chem. Phys., 2011, 13, 3960-3968] to real 
experimental conditions. These realistic suspensions include water dissociation ions and those gen- 
erated by atmospheric carbon dioxide contamination, in addition to the added counterions released 
by the particles to the solution. The electric potential at the particle surface will be calculated for 
different ion sizes and compared with classical Poisson-Boltzmann predictions for point-like ions, as 
a function of particle charge and volume fraction. The realistic predictions turn out to be essential 
to achieve a closer picture of real salt-free suspensions, and even more important when ionic size 
effects are incorporated to the electric double layer description. We think that both corrections have 
to be taken into account when developing new realistic electrokinetic models, and surely will help 
in the comparison with experiments for low-salt or realistic salt-free systems. 
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I. INTRODUCTION 

Many efforts have been devoted in the past, and still 
continue nowadays, with the aim of improving our knowl- 
edge of the electric double layer (EDL) surrounding a 
charged particle in a colloidal suspension [ll-Q- It is a 
well known fact that many non-equilibrium phenomena 
in this kind of systems are extremely sensitive to the 
specific properties of such EDL. For many years differ- 
ent electrokinetic models for colloidal suspensions have 
been derived. Most of them are based on the classi- 
cal Poisson-Boltzmann equation (PB), which is a mean- 
field theory with a reasonable success in representing the 
ionic concentration profiles at low to moderately charged 
interfaces in electrolyte solutions. However, one of its 
main drawbacks concerns the absence of size for the ions, 
which, for highly charged particles, yields to unphysi- 
cally high counterion concentration profiles near such in- 
terfaces. In addition, the PB treatment neglects ion-ion 
correlations, which simplifies the real scenario. Thus, the 
interaction on a particular ion is just represented by that 
in a mean-field, which is considered to be a poor descrip- 
tion of real interactions taking place inside the EDL by 
the microscopic models, in many cases [11,11]. 

On the other hand, most of the theoretical studies with 
colloidal suspensions corresponds to the dilute regime 
in particle concentration, in spite of it is the concen- 
trated one that deserves more interest due to its many 
industrial applications [7l-[l0j . The reason has to do with 
the larger complexity associated with the electrohydrody- 
namic particle-particle interactions in such concentrated 
systems, that are very difficult to manage theoretically. 
This has encouraged us to consider the latter particle- 



* E-mail address: carrique@uma.es 



particle interactions and the possibility of overlapping 
between adjacent double layers which will be unavoid- 
ably present with high particle concentrations. In many 
typical cases, the presence of an external salt added to the 
system gives rise to an effective screening effect on repul- 
sive electrostatic particle-particle interactions, depend- 
ing on the salt concentration, which are mainly respon- 
sible, for example, of the generation of colloidal crystals 
or glasses. Thus, it would be of worth to study systems 
with low screening regime for such interactions. 

Those systems are named salt-free because of the ab- 
sence of added external salt. The formation of colloidal 
crystals is easier in this kind of systems, even at suffi- 
ciently low particle volume fractions [HI, [l^l- Of course, 
these salt-free systems contain ions in solution, the so- 
called "added counterions" stemming from the particles 
as they get charged, that counterbalance their surface 
charge preserving the electroneutrality [l3l - [l7j . These 
salt-free suspensions have acquired a renovated interest 
in the last few years due primarily to the special phe- 
nomenology they show related to the colloid crystals. 

The authors have recently developed an EDL model 
of a spherical particle in a salt-free suspension and a 
robust and efficient mathematical treatment to numer- 
ically solved the PB equation, which for that case may 
become integrodifferential due to the coupling between 
particle charge and ionic countercharge in the solution. 
In addition, realistic conditions associated to the pres- 
ence of additional ions dissolved in the liquid medium, 
like those associated with water dissociation and possible 
atmospheric contamination, were included in the analy- 
sis [l3|. The aim was to improve the EDL model that 
could be used to develop nonequilibrium models in con- 
centrated systems, starting from the salt-free ones be- 
cause of their special theoretical interest. Thus, static 
and dynamic electrophoresis and complex conductivity 
and dielectric response models were developed by the au- 
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thors according to such EDL representation [l8l - [2l| . The 
first conclusion that can be drawn from these studies is 
the large magnitude of the realistic corrections for any of 
the latter properties. As a consequence, the neglecting 
of such corrections could lead to poor comparisons with 
experimental results. We can assure that the better the 
EDL representation of a real charged particle in a col- 
loidal suspension, the better predictions the macroscopic 
nonequilibrium models will be able to do. 

It is important to realize that the mean-field PB ap- 
proximation neglects ion-ion correlations, which may 
be responsible of some phenomena like overcharging or 
charge inversion that historically have not been predicted 
within this classical framework [22|. A first attempt to 
include some of these correlations concerns those linked 
to the ion excluded volume. The point-like ions of the 
PB approach are now considered as ions with finite size 
through a modified Poisson-Boltzmann equation (MPB) 
[23l - l39l | . It has been recently shown that these MPB pre- 
dictions of EDL equilibrium properties for dilute suspen- 
sions with a monovalent salt agree well with those ob- 
tained by Monte Carlo simulations or microscopic mod- 
els, although the agreement worsens as the valence of the 
counterions increases. Unfortunately, this latter MPB 
model showed to be unable to predict the overcharging 
phenomena whereas Monte Carlo simulations or other 
microscopic models succeeded @. 

Very recently, Lopez-Garcfa et al. [i^ have presented 
a mean-field MPB approach that includes finite ion size 
corrections with the additional assumption of a different 
distance of closest approach to the particle surface for 
counterions than for colons. They show that this model 
can predict charge inversion in the case of high elec- 
trolyte concentrations and counterion valence. We think 
that this is a very important result because it demon- 
strates that a mean-field theory can predict charge in- 
version, while microscopic models or Monte Carlo simu- 
lations achieve similar results by considering full ion-ion 
correlations. 

In this paper we will use an analogous MPB approach 
as the one previously derived for the pure salt-free case 
that included ionic size effects [4l|. It generalizes that 
of Borukhov [2^ for dilute salt-free systems to the con- 
centrated one, with the additional incorporation of an 
excluded region of a hydrated ion radius in contact with 
the particle. This latter issue has been shown by Aranda- 
Rascon et al. [s^ to provide a better representation of 
the solid-liquid interface and more reliable electrokinetic 
predictions. We will extend the finite ion size formalism 
to the case of realistic concentrated salt-free suspensions, 
already studied by the authors for point-like ions [l3| . 
As both corrections have considerably modified the stan- 
dard PB predictions, it is mandatory to know which are 
the predictions of a general mean-field approach includ- 
ing both, realistic salt-free corrections and finite ion size 
corrections. This is a difficult task because of the increas- 
ing numerical problems arisen when simultaneously take 
them into account in the resulting MPB equations de- 



rived for each case. It is also important to point out that 
the numerical instabilities progressively grow because of 
the iterative methods, unlike previous models for realistic 
salt-free systems where commonly no iterative methods 
were necessary. In the following sections we will present 
the theoretical model for this kind of systems. The re- 
sulting equations will be numerically solved and the equi- 
librium potential at the particle surface will be analyzed 
upon changing particle volume fraction, particle surface 
charge density, and size of the ions. In order to show 
the realm of the finite ion size effect in realistic salt-free 
suspensions, the results will be compared with MPB pre- 
dictions that do not take into account a finite distance 
of closest approach to the particle surface, and also with 
standard PB predictions for point-like ions. 



II. THEORY 

A. The cell model 

To account for the interactions between particles in 
concentrated suspensions, a cell model is used (bare 
Coulomb interactions among particles arc included in an 
average sense, but ions-induced interactions between par- 
ticles as well as ion- ion correlations, are ignored). For 
details about the cell model approach see the excellent 
review of Zholkovskij et al. [43]. This approach has 
been successfully used by the authors in the study of 
DC and AC electrokinetics and rheological properties of 
pure [l3 - [l7| and realistic [Tsl. [l8l - [2H salt-free suspensions 
with point-like ions. We have learned from those works 
how important the description of the EDL is for the non- 
equilibrium theoretical responses. 

Concerning the cell model. Figure [U each spherical 
particle of radius a is surrounded by a concentric shell of 
the liquid medium, having an outer radius h such that the 
particle/cell volume ratio in the cell is equal to the par- 
ticle volume fraction throughout the entire suspension. 




FIG. 1: Cell model in a realistic salt-free suspension includ- 
ing a distance of closest approach of the ions to the particle 
surface. 
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(2.1) 



The basic assumption of the cell model is that the 
macroscopic properties of a suspension can be obtained 
from appropriate averages of local properties in a unique 
cell. 



B. Finite size of the ions 

In a very recent study we addressed the EDL of a 
pure salt-free suspension taking into account the finite 
size of the counterions [IJ]- In this work we deal with 
realistic salt-free suspensions and, consequently, we will 
have more ionic species dissolved in the liquid medium, 
that will be coupled by appropriate chemical equilibrium 
mass-action equations. Our systems consist of aqueous 
suspensions deionized maximally without any electrolyte 
added during the preparation. Hence, in addition to the 
added counterions released by the particles to the solu- 
tion, we must also consider the and OH" ions from 
water dissociation. Moreover, if the suspension is open 
to the atmosphere, there will be other ions produced by 
the atmospheric COj contamination. All these new ionic 
species could be coincident, or not, with that of the added 
counterions. Here, we will describe the general case of N 
ionic species in the suspension. Details for each situation 
can be found in subsection III Dl 

Let us consider a spherical charged particle of radius 
a and surface charge density a immersed in a realistic 
salt-free medium with N ionic species including that of 
the added counterions, that counterbalance its surface 
charge. 

In our description, the axes of the spherical coordi- 
nate system (r, 9, tp) are fixed at the center of the par- 
ticle. In the absence of any external field, the particle is 
surrounded by a spherically symmetrical charge distribu- 
tion. 

Within a mean-field approximation, the total free en- 
ergy of the system, F = U — TS, can be written in 
terms of the equilibrium electric potential 'l'(r) and the 
ionic concentration ^^(r) of the different ionic species, 
I = 1, . . . , iV, of the suspension. The configurational in- 
ternal energy contribution U is 



U 



dr 



N 



N 



^Zieni(r)*(r) -^/iin,;(r) 



(2.2) 



The first term is the self-energy of the electric field, 
where eo is the vacuum permittivity, and is the relative 
permittivity of the suspending medium. The next term 
is the sum of the electrostatic energies of the different 



ionic species in the electrostatic mean field, and the last 
term couples the system to a bulk reservoir, where /i^ is 
the chemical potential of the ionic species i. 
The entropic contribution ~TS is 



dr 



^ n,;(r) / ni{r) 
^ In ' 



+ 1- 



(2.3) 



where is Boltzmann's constant, T is the absolute tem- 
perature, and TT,™"^ is the maximum possible ionic con- 
centration due to the excluded volume effect, defined as 
^max _ where V is the average volume occupied by 
an ion in the solution. For simplicity, we assume that all 
types of ions have the same size, and therefore will 
take the same value for all of them. The first term inside 
the integral is the sum of the entropy of the different ionic 
species, and the second one is the entropy of the solvent 
molecules. This last term accounts for the ion size effect 
that modifies the classical Poisson-Boltzmann equation 
and was proposed earlier by Borukhov et al. [2^ 

The variation of the free energy F = U — TS with 
respect to ^'(r) provides the Poisson equation 



JV 



V2vl/(r) = Vz.n,(r) 



(2.4) 



and the ionic concentration of the ionic species i is ob- 
tained by performing the variation of the free energy with 
respect to ^^(r), yielding 



"j(r) 



bi exp ^- 



e'I'(r 



N 

E- 



exp 



(2.5) 



where bi is an unknown coefficient that represents the 
ionic concentration of the species i where the electric po- 
tential is zero. 

Applying the spherical symmetry of the problem and 
combining Equations 12 .41 and 12 . 51 we obtain the modified 
Poisson-Boltzmann equation (MPB) for the equilibrium 
electric potential 



d2^'(r) 2d*(r) 



dr^ 



dr 



N 



exp 



zte^(r) 



N 



exp 



knT 



- 1 



(2.6) 



We need two boundary conditions to solve the MPB 
equation. The first one is 



d4'(r) 



dr 



= 



(2.7) 
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which derives from the electroneutrahty condition of the 
cell and the application of Gauss theorem to the outer 
surface of the cell. The second one is 







(2.8) 



that fixes the origin of the electric potential at r = 6. 

The MPB problem, Equations [111 [2Jl and [HI can be 
solved itcratively using the electroneutrahty condition of 
the cell and appropriate chemical reactions to find the 
unknown bi coefficients (see subsection IIIDI for details 
on bi calculation). This kind of problem can be solved in 
a better way by using dimensionless variables , which 
are defined as 



a kbI 



:b,: n" 



rewriting Equation 12.61 as 
5(2;) 



(2.9) 



_ d^^{x) 2 d*(.T) 



dx^ 



da; 



N 



1=1 



N 7 



(2.10) 



i=l 



where wc have defined the function g{x). If we differen- 
tiate it, after a little algebra, we find that 



N 



9\x) + ^ 



Zibi- 



i=l 



+ -^5^(^)^'(^)=0 (2.11) 

where the prime stands for differentiation with respect 
to X. In terms of the electric potential, Eauation l2.11l is 
rewritten as 



2 ~ 



*"'(a;) + -*"(x) -^'{x) + 'if'{x) *"(x) + -^'{x) 



N 



-2i*(2:) 



N 



*"(x) + -^'{x) 

X 



(2.12) 



Eauation l2.12l is a nonlinear third-order differential equa- 
tion that needs three boundary conditions to completely 



specify the solution. Two of them are provided by Equa- 
tions [5771 and [2^ which now read 



^'{h) = 0; *(/i) = 



(2.13) 



where h = {b/a) = 0"^/'^ is the dimensionless outer ra- 
dius of the cell. The third one specifies the electrical state 
of the particle, and can be obtained by applying Gauss 
theorem to the outer side of the particle surface r — a 



d«'(r) 



dr 



Its dimensionless form is 



*'(1) = -a 



(2.14) 



(2.15) 



Many theoretical and experimental studies have re- 
cently explored the use of either constant surface po- 
tential or constant surface charge boundary conditions. 
These conditions, although making the mathematical 
treatment simpler, represent only limiting or idealized 
cases. Many biological and artificial particles have their 
surface charge associated with some degree of dissocia- 
tion of functional groups which depend on the nearby 
environment. Constant surface potential and constant 
surface charge models would correspond, respectively, to 
the cases when the dissociation reactions of the func- 



tional groups are infinitely fast and infinitely slow [44 1. 
Also, some authors have derived an hybrid surface charge 
model to account for the electrical state of the parti- 
cles [13, which is a generalization of the conventional 
constant surface potential and constant surface-charged 
density models. Our model can be modified to include 
charge regulation mechanisms at the particle surface. 

It is very common in the literature to use the surface 
charge or the surface potential as a boundary condition at 
the particle surface when solving the equilibrium Poisson- 
Boltzmann equation, and both of them are valid. When 
it comes to concentrated suspensions, we prefer to use 
the particle surface charge boundary condition be- 
cause in many cases of interest the particle charge is a 
property that can be determined experimentally. More- 
over, the use of commercially available latex suspensions 
with fully dissociated surface electrical groups, have led 
us to choose the constant surface charge boundary condi- 
tion. In the near future our intention is to make electroki- 
netic predictions based on our model that can be checked 
against experimental data, and the latex suspensions are 
probably the most promising ones for that task. 

Besides, the surface potential depends on the choice of 
the potential origin, and in the case of concentrated sus- 
pensions there is not a standard criterium for this choice. 
As stated in Equation 12.81 wc have chosen it at the outer 
surface of the cell, r — b. 

If we consider point-like ions, = 00, Equation 

12.121 generates the expressions obtained by Ruiz-Reina 
and Carrique [l^ for realistic salt-free suspensions. Also, 
if we evaluate Equation l2 . 1 2l considcring that wc have just 
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one ionic species, the added counterions that counterbal- 
ance the particle surface charge, the equation reduces to 
the one obtained by the authors [4lj for pure salt-free 
suspensions including finite ion size effects. 

It can be easily demonstrated that solving the third- 
order problem is mathematically equivalent to finding the 
solution of the MPB problem, that is, any function sat- 
isfying Equations 12.121 12.131 and 12.151 also satisfies Equa- 
tions [2]6l [2?7l and 12.81 and vice versa. In both cases, we 
will use an appropriate iterative process for the calcula- 
tion of the bi coefficients in the resolution of the differ- 
ential equation. We set initial values for the different bi 
coeflScients, but the use of the third order problem with 
the addition of the boundary condition Equation 12.151 
improves the convergency giving rise to a better numeri- 
cal resolution. 

Equations 12.121 12.131 and 12.151 form a boundary value 
problem that can be numerically solved by using the 
MATL AB routine bvp4c [1^ , that computes the solution 
with a finite difference method by the three-stage Lobatto 
IIIA formula. This is a collocation method that provides 
a solution that is fourth order uniformly accurate at 
all the mesh points. The resulting mesh is non-uniformly 
spaced out and has been chosen to fulfill the admitted 
error tolerance (taken as 10^^ for all the calculations). 

Once we have found the electric potential ^(cc) with 
the iterative process for the calculation of the bi coeffi- 
cients, we can obtain the ionic concentration ni{r) for the 
different ionic species (see subsection III Dl for details). 



C. Excluded region in contact with the particle 

Following the work of Aranda-Rascon et al. [s^, we 
incorporate a distance of closest approach of the ions to 
the particle surface, resulting from their finite size. As we 
said before, for the sake of simplicity we have considered 
that all ions have the same size. We assume that ions 
cannot come closer to the surface of the particle than the 
chosen effective hydration ionic radius, _R, and, therefore, 
the ionic concentration will be zero in the region between 
the particle surface, r = a, and the spherical surface, 
r = a + defined by the ionic effective radius. 

The whole electric potential 5'(r) is now determined 
by combining Laplace's and MPB equations into the fol- 
lowing MPBL stepwise equation 



Equation 12.61 



a <r < a + R 
a + R<r <b 



(2.16) 



We must impose the continuity of the potential and of 
its first derivative at the surface r = a + R, in addition 
to boundary conditions. Equations 12 . 71 and 12.81 The con- 
tinuity of the first derivative comes from the continuity 
of the normal component of the electric displacement at 
that surface. Thus, in the region in contact with the par- 
ticle [a, a-|-i?], we are solving Laplace's equation, and, in 



the region [a + R, b], the MPB equation that we obtained 
previously in Equation 12.61 

As we have seen before, changing the system of second 
order differential equations. Equation 12.161 into one of 
third order, hugely simplify the resolution process 

*"'(a;) + f ^"(a;) - -^^'{x) =0 l<x <1 + S 

Equation [21^ l + S <x <h 

(2.17) 

where 6 = R/a and we have use dimensionless variables. 
The boundary conditions needed to completely close the 
problem are 

^'Lil) = -a ^'^(1 + 6) = 



*p(/i) = 



"^'pih) = 

^'^{1 + 5) ^ ^'p{l + S) 

(2.18) 



where subscript L refers to the region in which the poten- 
tial is calculated using Laplace's equation, and subscript 
P refers to the region in which we evaluate the MPB 
equation. 

The spherical solution of Laplace equation in the region 
free of charge is 



^l{x) = -+K 

X 



(2.19) 



where K is a constant. The potential difference in the 
Laplace layer becomes 



A^L = *l(1) -*i(l + 5) 



a5 



aS (2.20) 



where we have assumed that a ^ R. We solve now the 
MPB problem only in the [1 + S, h] region, instead of the 
[1, h] range that is used in the absence of a Laplace layer. 
Both solutions are very similar because 5 <^h and 5 <C 1 
and, consequently, we have 



^!p(l + S) w §(1) 



(2.21) 



where \t'(l) is the electric potential at the particle surface 
in the MPB case without any excluded region in contact 
with the particle. As a result, we expect that the surface 
potential will be shifted in accordance with Equation l2.201 
when we introduce the Laplace layer, in comparison with 
the MPB problem. The shift depends approximately lin- 
early on the Laplace region thickness 5 and, therefore, 
it remains roughly constant upon changing the volume 
fraction. 

The complete electric potential ^{x) is obtained nu- 
merically using Equations 12.171 and 12.181 including an 
iterative process for the calculation of the bi coefficients. 
Once the electric potential is found, the ionic concentra- 
tions riiir) for the different ionic species can be derived 
(see subsection III Dl for details). 
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D. Particularization of the MPBL equation for 
realistic salt-free concentrated suspensions 

In a previous work [4l| , the authors considered the case 
of pure salt-free concentrated suspensions inchiding ion 
size effects. In this subsection, we wih focus our treat- 
ment to rcahstic salt-free concentrated suspensions, with 
the consideration of ions from water dissociation and at- 
mospheric contamination, in addition to the added coun- 
tcrions. This will be carried out by using the generalized 
MPBL stepwise equation. Equation \2.1Ql or in a better 
way, its third order equivalent version given by Equation 
12.171 previously obtained. 



1. Added counterions and water dissociation 

Let us consider that, in addition to the added counte- 
rions stemming from the particle charging process, there 
are also H+ and OH" ions coming from water dissoci- 
ation in the liquid medium. This will always occur in 
aqueous suspensions. The equilibrium mass-action equa- 
tion for water dissociation, which we assume to hold at 
the outer surface of the cell, is 



[H+][OH-] = /w 6h+6oh- = 



(2.22) 



where the square brackets stand for the molar concen- 
tration, = 10""'^* mol^/L^ is the water dissociation 
constant at room temperature, 298.15 K, and A'w is a 
dimensionless quantity defined by 



i^w = ( —r- 1 

eoerksT 



(2.23) 



with Na the Avogadro constant. 

We can distinguish between two cases, (a) when the 
added counterions are or OH" ions, and (b) when 
they are of a different ionic species. The distinction is 
important because in the (a) case, the added counterions 
will enter in the equilibrium reaction equation for water 
dissociation, whereas in the (b) case, they do not. 

Case a. In this case we have two different ionic 
species. Evaluating Equation I2.10[ particularized just 
for H^ and OH" ions, a.t x ~ h we obtain 



use Equations 12.221 and 12.241 to find a new value b^]. , 

which will give us ^^^^(x) using Equations 12 . 1 71 and 12 . 181 
again. The numerical iterative process is repeated until 
the relative variation of the electric potential at the parti- 
cle surface is lower than a prescribed quantity. Although 
an iterative process has been used to obtain the solu- 
tion to Equation 12.171 with boundary conditions. Equa- 
tion 12.181 and Equations 12.221 and 12.241 this procedure 
is much better than the original and equivalent iterative 
problem defined by Equation l2.161 The improved conver- 
gency and superior numerical efficiency that are obtained 
when computing the third order problem lie in the facts 
that all of the intermediate solutions VP^"-* (x) of the iter- 
ative method have the correct slope ^'(1) = —a at the 
particle surface. This is not true if we use the original 
scheme because in that case the slope at the particle sur- 
face is not determined by any condition. 

Case b. In this case we have three different ionic 
species. Evaluating Equation l2.101 particularized for the 
added counterions, of valence Zc, and the ions H"*" and 
OH" , at x ~ h we obtain 

^"{h) = -zjbc - 2:h+^h+ - 2oH-^OH- (2-25) 

The added counterions counterbalance the overall 
charge on the particle surface 

-x^dx (2.26) 



1 ^ I 

Z J ^max y 

1=1 

whereas the number of H^ and OH" must be equal due 
to the electroneutrality of the cell. 

Using the relation between both 6h+ and 6oh- coef- 
ficients given by Equation 12.221 we can write Equations 
12.171 [2^251 and l2.26l in terms of two unknown coefficients, 
be and 6fj+ . The iterative process to obtain the solu- 
tion to Equation 12.171 with boundary conditions. Equa- 
tion [^TE and Equations and \TM is similar 
to the one described before, but in this case we have the 
two unknown coefficients, be and 6fj+ j to be determined 
iteratively. 



2. Added counterions, water dissociation, and atmospheric 
contamination 



-^h+Oh+ 



^OH-^OH- 



(2.24) 



Using the relation between both ^-^d 6oh- coef- 
ficients given by Equation 12.221 we can write Equations 
12.171 and 12.241 in terms of just one unknown coefficient, 
. The iterative process needed for the numerical res- 
olution of the third order version of the MPB problem 
remains as follows. We choose an initial guess for the 
&H+ coefficient, say b^l = or, in a more accurate way, 
the value obtained for point-like ions, and solve Equa- 
tion 12.171 with the boundary conditions given by Equa- 
tion We obtain the solution ^'^^^x), and then, we 



Let us consider now that, in addition to the added 
counterions and the H^ and OH " ions coming from wa- 
ter dissociation, there are also present ions stemming 
from the atmospheric COj contamination in the liquid 
medium. This will always occur in aqueous suspensions 
in contact with the atmosphere; the CO2 gas diffused 
into the suspension combines with water molecules to 
form carbonic acid HjCOg, and then, the following dis- 
sociation reactions take place 



H2CO3 
HCO^ 



H^ 
H^ 



HCO3 

cor 



(2.27) 
(2.28) 
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with equilibrium dissociation constants Ki = 4.47 • 10~^ 
mol/L and A'2 = 4.67 • 10^^^ mol/L at room tempera- 
ture, 298.15 K, respectively. The concentration of H2CO3 
molecules in water can be calculated from the solubility 
and the partial pressure of CO2 in standard air. For a 
temperature of 298.15 K and an atmospheric pressure 
of 101300 Pa, the concentration of carbonic acid is ap- 
proximately [H2CO3] = 1.08 • 10^^ mol/L, being its par- 
ticular value dependent on the local environmental con- 
ditions. The dimensionless dissociation constants and 
H2CO3 concentration are 



and HCOo , at x = h wc obtain 



~ lO^iV^e^^^ ~ 
1^1 = i-T^Ki] K2 



eoerksT 



eoerksT 



K2 

TVh.co, (2.29) 



where all the values are taken in S.I. units. The equilib- 
rium mass-action equations, which we assume to hold at 
the outer surface of the cell, are 



*"(/i) 



^OH-"OH- 



(2.33) 

The added counterions counterbalance the overall 
charge on the particle surface, as in Equation I2.26[ 
whereas the number of ions H"*", OH" and HCO3 must 
balance due to the electroneutrality of the cell. 

Using the relations between coefficients &h+i ^oh- ^-^d 
^HCO" given by Equations 12.221 and 12.301 we can write 
Equations 12.171 12.331 and 12.261 in terms of two unknown 
coefficients, be and 6h+ . The iterative process to obtain 
the solution to Equation 12. 171 with boundary conditions. 
Equation and Equations \m[ \TM 1^351 and [1211 

is similar to the case (b) for added counterions and water 
dissociation described before. 



III. RESULTS AND DISCUSSION 



[H+][HC03-] 
[H2CO3] 

[H+][C03 



[HCO3 



= K2 



^H+^CO., 
^HCOT 



= K- 



ki (2.30) 
(2.31) 



Hereafter, the second dissociation reaction. Equation 
12.281 will be neglected because the terms associated to 
the ion CO;^ that appears in Equation l2 . 1 71 are several or- 
ders of magnitude lower than those due to the ion HCO3" , 
in accordance with what the authors showed in a previous 
work [Tsj . 

Once more, we can distinguish between two cases, (a) 
when the added counterions are coincident with one of 
the ionic species in the system (H"*", OH" or HCO3 ) 
and (b) when they are of a different ionic species. In 
the (a) case, the added counterions will enter in one of 
the equilibrium dissociation equations, whereas in the (b) 
case, they do not. 

Case a. In this case wc have three different ionic 
species. Evaluating Equation l2.101 particularized for H^, 
OH" and HCO3 ions, &i x — h we obtain 



'^OH-"OH- 



(2.32) 



Using the relations between coefficients &h+ i ^oh^ 
6jjQQ- given by Equations 12.221 and 12.301 wc can write 
Equations 12.171 and 12.321 in terms of just one unknown 
coefficient, 6h+- The iterative process to obtain the solu- 
tion to Equation 12.171 with boundary conditions, Equa- 
tionUm and Equations l^T^IOUl and E3^ is similar to 
the case (a) for added counterions and water dissociation 
described before. 

Case b. In this case we have four different ionic 
species. Evaluating Equation l2.101 particularized for the 
added counterions, of valence Zc, and the ions H^, OH" 



For all the calculations, the temperature T has been 
taken equal to 298.15 K and the relative electric permit- 
tivity of the suspending liquid = 78.55, which coin- 
cides with that of the dcioniscd water. Also, the valence 
of the added counterions Zc has been chosen equal to -|-1, 
when they are of a different ionic species as that of the 
ions stemming from water dissociation or atmospheric 
CO2 contamination, and the particle radius a = 100 nm. 
Other values for Zc could have been chosen. The model 
for point-like ions is able to work with any value of Zc in 
the Poisson-Boltzmann equation, but we think that the 
predictions of this model will be less accurate in the case 
of multivalent counterions, since it is based on a mean- 
field approach that does not consider ion-ion correlations, 
which are increasingly more important as the ion valence 
grows. Nevertheless, when we take into account the fi- 
nite size of the ions, the main objective of this work, we 
include correlations associated with the ionic excluded 
volume, solving partially this problem because we are 
still not considering electrostatic ion-ion correlations. 

For the sake of simplicity, we assume that the average 
volume occupied by an ion is ^ = (2i?)^, being 2R the 
chosen effective ionic diameter. With this consideration, 
the maximum possible ionic concentration due to the ex- 
cluded volume effect is n"^"^^ = {2R)~^. This corresponds 
to a cubic package (52% packing). In molar concentra- 
tions, the values used in the calculations, = 22, 4 
and 1.7 M, correspond approximately to effective ionic 
diameters of 2R = 0.425, 0.75 and 1 nm, respectively. 
These are typical hydratcd ionic diameters (47j . 

We will discuss the results obtained from three dif- 
ferent models, the classical Poisson-Boltzmann equation 
(PB), the modified Poisson-Boltzmann equation by the 
finite ion size effect (MPB), and the MPB equation in- 
cluding also a distance of closest approach of the ions 
to the surface of the charged particle (MPBL). Besides, 
we will consider different realistic salt-free concentrated 
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FIG. 2: Dimensionless surface electric potential as a function 
of particle volume fraction, considering (blue lines, MPBL) 
or not (red lines, MPB) the excluded region in contact with 
the particle. Black lines show the results of the PB equation. 
Solid lines: only added counterions H"'' (SF). Dashed lines: 
added counterions H"'" with water dissociation ions (WDC). 
Dotted lines: added counterions different than with Zc — 1 
and water dissociation ions (WDNC). 



suspensions: suspensions with added counterions coinci- 
dent, WDC (read it as Water Dissociation Coincident), or 
not, WDNC (Water Dissociation Non-Coincident), with 
the ions that stem from water dissociation, i.e., H'*' or 
OH" ; and suspensions with added counterions coinci- 
dent, WDACC (read it as Water Dissociation Atmo- 
spheric Contamination Coincident), or not, WDACNC 
(Water Dissociation Atmospheric Contamination Non- 
Coincident), with the ions that stem from water disso- 
ciation or atmospheric contamination, i.e., H^, OH", 
HCO3 . We will compare them with the results of pure 
salt-free suspensions with only added counterions (SF), 
obtained by the authors in a previous work [4l| . 

Fig. [5] shows the dimensionless equilibrium electric po- 
tential at the surface of the particle for a wide range of 
particle volume fractions. We display in black, red and 
blue hues the predictions of the PB, MPB and MPBL 
equations, respectively. Solid lines account for the results 
of pure salt-free suspensions with only added counteri- 
ons H"*" (SF). Dashed hues stand for realistic WDC sus- 
pensions with H^ ions as added counterions, and dotted 
lines for realistic WDNC suspensions with added counte- 
rions different than H"*" with Zc = 1. The particle surface 
charge density have been chosen equal to —40 /iC/cm^, 
and the maximum possible ionic concentration due to the 
excluded volume effect = 22 M. 

We can observe in Fig. [5] that for high particle con- 
centration, the results for realistic suspensions with wa- 
ter dissociation ions are equal to those for pure salt-free 
suspensions, because the added counterions completely 
mask the influence of water dissociation ions, due to its 
larger concentration. The situation is very different when 



we approximate the dilute limit. For WDC suspensions 
with particle volume fraction (j) < 10~^ the surface po- 
tential becomes approximately constant, in contrast with 
the growth noticed in a pure salt- free suspension (SF). 
This behaviour is due to the different sources of H"*" coun- 
terions: the added counterions released by the particles 
and those stemming from water dissociation. The first 
ones dominate in the high region, whereas the second 
ones do it at low particle concentration. The diminution 
of the surface potential in comparison with the pure salt- 
free case can be explained by the increase of the coun- 
terion concentration inside the cell. For WDNC suspen- 
sions (dotted Hues), there is an additional decrease of the 
surface potential and a wider influence of water dissocia- 
tion to higher volume fraction values in comparison with 
the case of coincident counterions. WDC (dashed lines). 
This is due to the fact that now the added counterions do 
not participate in the water dissociation reaction. There- 
fore, the total number of counterions is larger than in 
the WDC case, causing a better screening of the particle 
charge and, consequently, diminishing the electric poten- 
tial at the particle surface. 

Fig. [2] also shows that the inclusion of the finite ion size 
effect (MPB calculations) always rises the surface electric 
potential in comparison with the PB predictions. The 
reason relies on the limitation of the ionic concentration 
in the neighborhood of the particle, which seriously di- 
minishes the screening of the particle charge, and conse- 
quently, leading to an increment of the surface potential. 
An additional significative increase of the surface poten- 
tial is obtained when the distance of closest approach 
of the ions to the particle surface is taken into account 
(MPBL model). As expected, the existence of a Laplace 
region free of ions also penalizes the screening of the par- 
ticle charge. We always find this behaviour whatever the 
cases studied: SF, WDC or WDNC suspensions. As we 
mentioned before, there is a potential shift between the 
red and blue lines which is independent of the volume 
fraction and also agrees numerically well with Equation 
12.201 This explain the parallelism observed between these 
curves. 

In the MPB case without a Laplace layer, and for suf- 
ficiently highly charged particles, a region of constant 
charge density develops very close to the particle surface 
with a thickness that is approximately independent of 
the volume fraction. This is in contrast with the classi- 
cal PB problem, where the charge density at the particle 
surface is unbounded and can reach unphysical values. It 
is clear that we can apply a similar explanation as before 
for the parallelism that also exists between the red and 
black curves. 

Fig. Oshows the dimensionless equilibrium electric po- 
tential at the surface of the particle for a wide range of 
particle volume fractions. We display in black, red and 
blue lines the predictions of the PB, MPB and MPBL 
equations, respectively. Solid lines account for the re- 
sults of pure salt-free suspensions with only added coun- 
terions H^ (SF). Dashed lines stand for realistic WDACC 
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FIG. 3: Dimensionless surface electric potential as a function 
of particle volume fraction, considering (blue lines, MPBL) 
or not (red lines, MPB) the excluded region in contact with 
the particle. Black lines show the results of the PB equation. 
Solid lines: only added counterions H"*" (SF). Dashed lines: 
added counterions with water dissociation and atmo- 
spheric contamination ions (WDACC). Dotted lines; added 
counterions different than H"*" with Zc = 1 and water dissoci- 
ation and atmospheric contamination ions (WDACNC). 



FIG. 4: Dimensionless surface electric potential for different 
values of particle surface charge density (a) and particle vol- 
ume fraction (b), considering (blue lines, MPBL) or not (red 
lines, MPB) the excluded region in contact with the parti- 
cle. Black lines show the results of the PB equation. Solid 
lines: only added counterions H'*' (SF). Dashed lines: added 
counterions with water dissociation ions (WDC). Dotted 
lines: added counterions H"*" with water dissociation and at- 
mospheric contamination ions (WDACC). 



suspensions with H ions as added counterions, and dot- 
ted lines for realistic WDACNC suspensions with added 
counterions different than with = 1- In the case 
of WDACC (dashed lines), we again observe that for 
highly concentrated suspensions the results coincide with 
those of pure salt- free predictions (SF), due to the added 
counterions dominance over water dissociation and atmo- 
spheric contamination ions. In contrast with that shown 
in Fig. [5] for water dissociation ions, when we also con- 
sider atmospheric CO2 contamination, the plateau in the 
surface potential now extends from the very dilute limit 
to ^ = 10~^. In the present case, the counterions in- 
side of the cell arise from three different mechanisms: the 
charging process of the colloidal particle, the water dis- 
sociation equilibrium, and the dissociated protons from 
atmospheric carbonic acid. Consequently, there is a great 
increase of the concentration of counterions that accounts 
for the marked reduction of the surface potential in the 
low volume fraction region. 

For WDACNC suspensions (dotted lines). Fig. H dis- 
plays an additional decrease of the surface potential and 
an extended influence of the atmospheric CO2 contami- 
nation for larger volume fractions in comparison with the 
case of WDACC suspensions. The explanation is based 
again on the fact that the non coincident added counte- 
rions do not participate in the water and carbonic acid 
dissociation reactions, yielding a large number of counte- 
rions in solution. The resulting screening of the particle 
charge is enhanced in comparison with the WDACC case, 
and therefore, the surface potential decreases. 

Regarding the finite ion size effect. Fig. [3] shows that 



the surface potential always increases for MPB calcula- 
tions in comparison with the PB predictions. This is 
again due to the limitation of the ionic concentration in 
the neighborhood of the particle surface, which provokes 
a diminution of the screening of the particle charge, and 
consequently, raising the surface potential. The surface 
potential increases in a significative way when we con- 
sider a distance of closest approach of the ions to the par- 
ticle surface (MPBL model). The reason lies on the re- 
duction of the particle charge screening as a consequence 
of the Laplace region free of ions next to the surface. We 
again find this behaviour irrespective of the cases studied: 
SF, WDACC or WDACNC suspensions. 

Fig. 2] displays the dimensionless equilibrium sur- 
face electric potential for a wide range of particle sur- 
face charge densities. Fig. |4^, and particle volume frac- 
tions. Fig. lib- We repeat this study for different salt- 
free suspensions, SF, WDC and WDACC, with added 
counterions coincident with H"*", and for the three PB, 
MPB and MPBL models. We take the maximum possi- 
ble ionic concentration due to the excluded volume effect 
as 71™°^ = 22 M, corresponding to a hydrated hydro- 
nium ion in solution. We observe in Fig. 2^ that the 
diminution of the surface potential due to realistic con- 
siderations (dashed or dotted lines against solid lines) 
is practically independent of the particle surface charge 
density, for both point-like and finite size ions, even when 
considering the excluded region in contact with the par- 
ticle. The results of Fig. confirm those from Figs. [2] 
and[21 and we can clearly see how the plateaus in the sur- 
face potential extend to larger volume fractions when we 
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FIG. 5: Dimensionless surface electric potential for different 
values of particle surface charge density (a) and particle vol- 
ume fraction (b) for different ion sizes, considering (dashed 
lines, MPBL) or not (solid lines, MPB) the excluded region 
in contact with the particle. Black lines show the results of 
the PB equation. In all cases we consider salt-free suspen- 
sions with added counterions different than H"*" with Zc = 1 
and water dissociation and atmospheric contamination ions 
(WDACNC). 



consider WD AC C (dotted lines) instead ofWDC (dashed 
lines) suspensions. 

Fig. [5] displays the dimensionless equilibrium surface 
electric potential for a wide range of particle surface 
charge densities. Fig. [5^, and particle volume fractions. 
Fig. We repeat this study for different ion sizes. In 
all cases we examine realistic salt-free suspensions with 
added counterions different than H"'' with Zc ~ 1 and 
water dissociation and atmospheric contamination ions 
(WDACNC). We find that the surface electric poten- 
tial increases with the particle charge density, Fig. [S^. 
However, in some cases there is a different behaviour in 
comparison with the point-like case. Initially, a fast and 
roughly increase of the surface potential with the sur- 
face charge density is observed, which is followed by a 
much slower growth at higher surface charge densities for 
the PB case, or when the size of the ions is very small. 
This phenomenon is related to the classical countcrion 
condensation effect: for high surface charges a layer of 
counterions develops very close to the particle surface 
0. When the ion size is taken into account (MPB), 
we limit the appearance of the classical condensation ef- 
fect because when the surface charge is increased, the 
additional counterions join the condensate enlarging it. 
If a region of closest approach of the ions to the par- 
ticle surface is also considered (MPBL), the mechanism 
is the same, but now the additional counterions are lo- 
cated in farther positions from the particle surface. This 
explains the additional increase of the surface potential 
observed for large ion sizes and high surface charges den- 
sities. These profiles are similar to those found by the 



FIG. 6: Dimensionless surface potential against the surface 
charge density for different particle volume fraction values. 
Solid lines stand for the results of the PB equation. Dashed 
lines show the results of the MPBL model. In all cases 
we consider salt-free suspensions with added counterions H"*" 
and water dissociation and atmospheric contamination ions 
(WDACC). 



authors for pure salt- free suspensions [41| , but the values 
of the surface potential are lower for WDACNC suspen- 
sions than for SF suspensions, as Fig. [3] displayed. 

On the other hand, the surface electric potential de- 
creases when the particle volume fraction increases, irre- 
spective of the cases studied: PB, MPB or MPBL, Fig. 
[SJd. When the particle concentration raises, the available 
space for the ions inside the cell decreases and, conse- 
quently, the screening of the particle charge largely aug- 
ments, thus reducing the value of the surface potential. 

Figs. [6] and [7] expand the results of Fig. |4] for real- 
istic salt-free suspensions with added counterions coin- 
cident with and water dissociation and atmospheric 
contamination ions (WDACC). In them, we compare the 
results of the PB equation (solid lines) with those of the 
MPBL model (dashed lines) for a given ion size (typi- 
cal of an hydronium ion in solution). Fig. [5] presents 
the dimensionless equilibrium surface electric potential 
at a wide range of particle surface charge densities. The 
different coloured lines correspond to different particle 
volume fraction. The most remarkable fact shown in 
Fig. [6] is the large influence of the finite ion size effect 
even for moderately low particle surface charge densities. 
While for PB predictions (solid lines) the surface poten- 
tial hardly increases with surface charge for moderate to 
high surface charges at each volume fraction, the MPBL 
results (dashed lines) display an outstanding growth for 
the same conditions, associated with the lower charge 
screening ability of the counterions because of their finite 
size. This fact will surely have important consequences 
on the electrokinctic properties of such particles in con- 
centrated realistic salt-free suspensions, as has already 
been shown for dilute suspensions in electrolyte solutions 
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FIG. 7: Dimensionless surface potential against the parti- 
cle volume fraction for different surface charge density val- 
ues. Solid lines stand for the results of the PB equation. 
Dashed lines show the results of the MPBL model. In all cases 
we consider salt-free suspensions with added counterions 
and water dissociation and atmospheric contamination ions 
(WDACC). 



by Aranda-Rascon et al. [36|, |37| . 

Finally, Fig. [7] shows the dimensionless equilibrium 
electric potential at the particle surface against the par- 
ticle volume fraction. The different coloured lines corre- 
spond to different negative particle surface charge densi- 
ties. As previously stated, the particle surface potential 
shows an initial plateau in the dilute region followed by a 
monotonous decrease with volume fraction at fixed par- 
ticle charge density. For the ionic size chosen, the larger 
the surface charge, the larger the relative increase of the 
surface potential at every volume fraction, due to the en- 
larging of the counterion condensate in the neighborhood 
of the particle surface. 

From the results, we think that it is clear that the in- 
fiuence of the finite ion size effect on the EDL description 
cannot be neglected for many typical particle charges and 
volume fractions, either in a pure or in a realistic salt-free 
suspension. 



electric double layer of a spherical particle in a realis- 
tic salt-free concentrated suspension, in an attempt to 
get us closer to real systems. These realistic suspensions 
include water dissociation ions and those generated by 
atmospheric carbon dioxide contamination, in addition 
to the added counterions released by the particles to the 
solution. The resulting model is based on a mean-field 
approach that has reasonable succeeded in modeling clcc- 
trokinetic and rhcological properties of concentrated sus- 
pensions. 

We have used a cell model approach to account 
for particle-particle interactions, and derived modified 
Poisson-Boltzmann equations (MPB and MPBL) which 
include such ion size effects. The theoretical procedure 
has followed that by Borukhovjl^ with the additional 
inclusion of an excluded region of closest approach of the 
ions to the particle surface [s^ H^]- The results have 
shown that the finite ion size effect (MPB equation) has 
to be taken into account for moderate to high particle 
charges at every particle volume fraction, and even more 
if a distance of closest approach of the ions to the particle 
surface is considered (MPBL model), irrespective of the 
realistic models (WDC, WDNC, WDACC or WDACNC) 
used. 

The equilibrium model presented in this paper will 
be used to develop nonequilibrium models of the re- 
sponse of a realistic salt-free concentrated suspension to 
external electric fields. Experimental results concerning 
the DC electrophoretic mobility, dynamic electrophoretic 
mobility, electrical conductivity and dielectric response, 
should be compared with the predictions of these future 
models to test them. To carry out such comparisons 
highly charged particles like those of some sulfonated 
polystyrene latexes might be used. These theoretical and 
experimental tasks will be addressed by the authors in 
the near future. 
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